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ABSTRACT 

Recent studies have demonstrated that detailed monitoring of gravitational microlens- 
ing events can reveal the presence of planets orbiting the microlensed source stars. 
With the potential of probing planets in the Galactic Bulge and Magellanic Clouds, 
such detections greatly increase the volume over which planets can be found. This 
paper expands on these original studies by considering the effect of planetary phase 
on the form of the resultant microlensing light curve. It is found that crescent-like 
sources can undergo substantially more magnification than a uniformly illuminated 
disk, the model typically employed in studying such planets. In fact, such a circularly 
symmetric model is found to suffer a minimal degree of magnification when compared 
to the crescent models. The degree of magnification is also a strong function of the 
planet's orientation with respect to the microlensing caustic. The form of the magni- 
fication variability is also strongly dependent on the planetary phase and from which 
direction it is swept by the caustic, providing further clues to the geometry of the 
planetary system. As the amount of light reflected from a planet also depends on its 
phase, the detection of extreme crescent-like planets requires the advent of 30-m class 
telescopes, while light curves of planets at more moderate phases can be determined 
with today's 10-m telescopes. 

Key vifords: Gravitational Lensing - Planetary Systems 



1 INTRODUCTION 

Recently, there has been a rapid growth in the number of 
extra-solar planets identified (Vogt et al. 2000; Marcy et al. 
2000). Driven by technological advances, most have been 
found by the identification of a 'Doppler reflex', although 
more novel techniques, such as planetary occultation (Char- 
bonneau 2000), are also proving fruitful. Collier Cameron 
et al. (1999) claim to have directly detected the starlight 
reflected from the planet orbiting r Bootis (although see 
Charbonneau et al. 1999). While it now appears that this r 
Bootis measure has not been confirmed in follow-up observa- 
tions, the identification of planets due to reflected starlight 
has a firm theoretical foundation (Sudarsky et al. 2000). 

All these techniques focus on stars in the vicinity of 
the sun, over a region of several tens of parsecs. At larger 
distances, however, planets can be detected orbiting the 
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compact objects responsible for gravitationally microlens- 
ing stars in the Galactic Bulge and Magellanic Clouds via 
their perturbative influence on the microlensing light curve 
(Wambsganss 1997). More recently, Graff & Gaudi (2000) 
and Lewis & Ibata (2000) (hereafter LI2000) have considered 
instead the detection of planets orbiting the microlensed 
sources. Using the models for 'hot Jupiter' type planets 
as a basis for their studies, these groups demonstrated 
that the small fraction of light reflected from the planet 
{up to X 10~*I/,) can be significantly magnified, resulting 
in an observable ~ 1 per cent deviation in the microlensing 
light curve. LI2000 also demonstrated how the monitoring 
of polarization through a microlensing event can probe the 
composition of the planetary atmosphere. 

In studying the microlensing light curve of planets or- 
biting the source stars, these original studies assumed that 
the planet undergoing microlensing can be represented by 
a circularly symmetric source. This situation only occurs 
at opposition, and at any other point in its orbit the planet 
will display phases akin to the moon and the inferior planets. 
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Figure 1. The planetary phases and orientations considered in 
this paper. The cross-lines and arrow at the top of the figure in- 
dicate the caustic and its direction of motion, with the cross-linos 
denoting the high magnification region of the caustic. Through- 
out the text the models are referred to as {i,6), where i is the 
model label (a,b or c) and 9 is the orientation. 



This paper considers the influence of these deviations from 
circular symmetry on the form of the microlensing event 
with the goal of examining their influence on the detectabil- 
ity of planets via gravitational microlensing. Section ^ de- 
scribes the numerical approach taken in this paper, and out- 
lines both the microlensing model and the choice of plane- 
tary phases. Section |^ describes the results of this study, 
including an examination of how these results influence the 
detectability of the planet, while Section ^ presents the con- 
clusions. 



2 METHOD 

2.1 Microlensing Model 

Several reviews on the physical aspects of microlensing in 
the Galactic Halo have recently appeared (e.g. Paczyhski 
1986) and this background material will not be reproduced 
here. As with the studies of Graff & Gaudi (2000) and LI2000 
this paper assumes that the planetary system is microlensed 
by a binary system. Such systems can result in extended re- 
gions of strong magniflcation in the vicinity of 'fold caustics'. 
Due to their extended nature, such caustics have proved use- 
ful probes of not only the surfaces of stars (Sasselov 1998; 
Gaudi & Gould 1999), including the detection of stellar spots 
(Han et al. 2000), but also the structure at the heart of 
quasars (Fluke et al. 1999). As with Graff & Gaudi (2000) 
and LI2000, it is assumed that the planet is swept by a fold 
caustic during a binary microlensing event. 

The magnification induced as a point source is crossed 
by a fold caustic is given by 



zH{x - Xc) + Mo 



(1) 



Here, x—Xc is the distance from the source to the caustic and 
H{x) is the Heaviside step function (Schneider et al 1992). 
The flux factor k represents the 'strength' of the caustic 
and Ho is the ampliflcation of all other images of the source 



far from the critical curve of the caustic (Chang 1984; Witt 
1990; Wook et al. 1998) 

In general, analytically calculating the light curve of an 
arbitrary source as it is microlensed by a fold caustic is non- 
trivial. Numerically, however, the situation is much more 
tractable. For this study, the magnification distribution due 
to the fold caustic is defined on a grid of pixels. The grid 
is chosen such that the pixel size is much smaller than the 
scale of the magnification map or any planetary structure. 
In considering pixels. Equation ^ can written (Lewis & Belle 
1998) 



Aa; 



- Aa; - ^/xl) + /xo 



(2) 



Here, Ax is the extent of the pixel, and Xk is the distance 
of the pixel to the caustic. In the following study we assume 
/lo = 1, although we will concentrate on the magnification 
due to the caustic. 



2.2 Planetary Phase 

Figure |l| summarizes the planetary phases and orientations 
employed in this study. The lower panel presents the models 
for the phase; running from (a) to (c) these go from a quar- 
ter moon configuration to an extreme crescent. In the upper 
part of the Figure, the relative orientation of the planetary 
phase to the incoming caustic (denoted by the cross-lined 
region and arrow) is presented. Each model, therefore, is 
labelled as {i,9), where i represents the phase (a, b or c) 
and 9 is the orientation. The crescent is assumed to be uni- 
formly illuminated, neglecting the effects of limb darkening. 
In studies of stellar microlensing, limb darkening has been 
found to produce a variations of several percent during the 
peak of a microlensing event (e.g. Yock 2000) and can be 
considered a perturbation to the results presented here. 

In comparing to previous work it should be noted that 
due to symmetry considerations, the magnification of model 
(a, 0) - henceforth referred to as the "uniform, circular 
source" - is identical to that of the circular source of the 
same radius. Note also that, to allow a comparison, in the 
following the "light curves" present the magnification over 
time, equivalent to assuming that the flux received from each 
model of the planet is the same, irrespective of the phase of 
the planet. The effect of planetary phase on the observed 
ffux, and hence detectability of t he planetary microlensing 
event, is addressed in Section 



3.3 



2.3 Numerical Approach 

In the numerical approach adopted, each source was gen- 
erated as depicted by models (a) , (b) and (c) , in Figure |l| 
As noted previously, the magniflcation pattern was laid out 
on a grid of pixels. Similarly, the light distribution over the 
planets, as illustrated in Figure 1, was also defined on a 
grid, taking account of the orientation with respect to the 
incoming caustic. The grid scale was fine enough to ensure 
insignificant numerical errors, with a comparison to analytic 
tests. In generating the lensed view of the source, the plan- 
etary profiled was convolved with the magnification map, 
providing the light curve as the source is swept by t he c aus- 
tic. Note that all light curves discussed in Sections 3.1 and 
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Figure 2. The light curves for model (a). The different line styles, described in the key, refer to the fiducial orientations displayed in 
Figure 1. The axes represent the magnification and time scale for the microlensing of a hot Jupiter-like system located in the Galactic 
bulge, allowing comparison with LI2000. 



3.2 are normalized such that the unlensed view of the source 
has unit flux. 



3 RESULTS 

3.1 Quarter Moon Illumination: Model (a) 

Figure ^ presents the light curves for model (a) at the flve 
fiducial orientations described in Section 2.2, The axes are 



for the model of LI2000 who considered a planet of radius 
1.8x Jupiter located in the Galactic Bulge, a value typical 
of hot Jupiter systems. The light curves presented in this 
paper can related to those of differing source size, caustic 
strength and relative velocities through simple scaling re- 
lations. Firstly, both the distances and caustic strength in 
Equations |l] and ^ are naturally expressed in units of Ein- 
stein radii. Given the relative velocities, the crossing time of 
this scale length is given by (e.g. LI2000); 



tE 



140 

^^100 



M Dos 
Mq 8kpc 



(1 — d)d days 



(3) 



where M is the mass of the lens (typically taken to be IMq) 
and Dioo is the perpendicular velocity of the lensing objects 
in units of lOOkms"^. Dos is the distance to the source in 
kpc, while d is the ratio of the distance to the lens to the 
distance to the source. The magnification can be scaled using 
the results of Chang (1984) who demonstrated that when a 
source of radius Rs is swept by a caustic of strength k, it 
will suffer a maximum magnification of 



^max — 



(4) 



where / is the 'form factor' which accounts for the specific 
shape and brightness profile of the source. For a uniformly 
illuminated circular disk, / = 1.39. 

It is immediately apparent that rotating the source with 
respect to the caustic radically changes the form of the mi- 
crolensing light curves. At (a, 7r/2), the light curve is seen 
to possess a rapid rise, sharply turning over into a rapid de- 
cline. Similarly, (a, — 7r/2) also displays a rapid rise, but a 
more gentle turn over into a decline. The intermediate cases 
show similar trends between the two extreme orientations. 

An examination of the light curves in Figure ^ reveals 
that the peak magnification in each light curve also depends 
on the planetary orientation. In investigating this further, 
one hundred additional light curves at orientations between 
— 7r/2 and 7r/2 were generated. Figure |^ presents the peak 
of the magnification of these light curves versus the orienta- 
tion of the planet with respect to the caustic (note that the 
background magnification, fio, has been subtracted from the 
light curves in this analysis). The curve is normalized with 
respect to the peak magnification of model (a,0), the uni- 
form, circular source (equivalent to a circular source of the 
same radius). Due to the scaling relations outlined previ- 
ously, this curve is insensitive to the specific values of the 
caustic strength and planetary radius. 

Interestingly, model (a) suffers significantly more mag- 
nification at orientations of ±7r/2 when compared to an ori- 
entation of 0; in the case of 7r/2 this is ~ 1.6x that of the 
uniform, circular source. This result clearly demonstrates 
that such a model undergoes a minimum magnification when 
compared to model (a) at virtually all other orientations. 
The degeneracy between model (a, 0) and a circularly sym- 
metric source of the same radius allows further comparison 



© 0000 RAS, MNRAS 000, 000-000 



4 C. E. Ashton & G. F. Lewis 




-0.4 -0.2 0.2 0.4 



Orientation / tt 

Figure 3. The dependence of peak magnification on tlie orien- 
tation of the planet for model (a). The peak magnifications are 
normalized to the light curve of the uniform, circular source; this 
will scale to previous investigations on circular sources, and allow 
comparisons. 
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Figure 5. The dependence of peak magnification on the orien- 
tation of the planet, normalized by that of the uniform, circular 
source. Model (b) is denoted by a solid line, and model (c) by a 
dotted line. 



with previous results. Hence, using Equation ^ the form fac- 
tor / for model (a) at the various orientations presented in 
this paper can be simply calculated by scaling by the values 
presented in Figure^ [using Figure]^ for models (b) and (c)]. 



3.2 Crescent Illumination: Model (b) & (c) 

Figure ^ presents the light curves for models (b) and (c). In 
comparing to Figure ^ it can be seen that the peak magni- 
fication of all models at 6' ~ is virtually the same. Both 
sets of light curves possess the same generic features seen in 
model (a), although as the planet becomes more crescent- 
like, the features of the light curve become sharper. This 
is strikingly apparent for model (c) (lower panel), where 
the light curve displays an almost vertical rise in (c, — 7r/4) 
and (c, — 7r/2). Correspondingly, (c, 7r/4) and (c, 7r/2) pos- 
sess substantial peaks in their light curves, although these 
are accompanied by more gentle rises. 

Again, the peak magnification is a function of the plan- 
etary orientation with respect to the caustic. As with model 
(a), a larger sample of light curves for models (b) and (c) 
were generated to investigate this dependence. The peak 
magnifications for these samples are presented in Figure ^ 
Interestingly, while it is apparent that light curves for each 
model at 6^ = are different, they possess very similar values 
of peak magnification. These have been normalized with the 
peak magnification of the uniform, circular source. As with 
model (a) , the light curves at 6' ~ possess the lowest value 
of peak magnification, with higher values occurring at ±7r/2, 
although in model (b) the extremes of the peak magnifica- 
tion are at ~ 1.6x the minimum value. Model (c) possesses 
a similar value at 7r/2, while at —it/2 the peak magnification 
climbs to ~ 2.4x the minimum value. 



3.3 Detectability 

The preceding sections have discussed the variation in mag- 
nification as a planet is swept by a gravitational microlensing 
caustic. While this eases comparison with previous models, 
being dependent only upon the size and shape of the source, 
the observed property during a microlensing event is fiux. 
For a specific planetary system, the fraction of stellar light 
refiected from a planet is dependent upon its orbital phase. 
As demonstrated in the previous sections, the details of the 
microlensing magnification are similarly dependent upon the 
planetary orbital phase, and it is the combination of the re- 
fiected planetary fiux and magnification that will determine 
detectability of the microlensing event. 

With respect to a circular source, models (a), (b) and 
(c) reflect 50%, 29% and 5% of the stellar light respectively, 
assuming isotropic scattering. Figure ^presents the observed 
fiux for a specific planet at several different phases, including 
gibbous phases (denoted in the inset box). Four orientations 
with respect to the sweeping caustic are portrayed, and all 
light curves are normalised to the fiux of a microlensed circu- 
lar source, which peaks at 100%. The same feature is seen in 
all light curves, namely that the observed fiux drops as the 
planetary phase changes from full to crescent, even though 
these latter phases undergo stronger magnification. 

Figure ^ presents simulated observations of light curves 
given in Figures 2 and 4, assuming the host star is V=18, 
with an extinction of 1.25mags (following the model of 
LI2000). The three models, (a), (b) and (c) are presented 
vertically, with each column representing the differing ori- 
entations. For each model, two sets of simulated observations 
are given, assuming a 10-m class (upper panels) and 30-m 
class (lower panels), using the planetary parameters used by 
LI2000. The ordinate on each panel represents the fractional 
light curve deviation due to the presence of the planet, and 
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Figure 4. As Figure 2, but for models (b) and (c). The orientations in the key refer to the fiducial models denoted in Figure 1. In 
comparing to Figure ^ it should be noted that the peak magnification of all models at 6 ^ is virtually the same. As with Figure 2, 
considering the microlensing of a hot Jupiter-like planet in the Galactic Bulge. 



should be compared to Figure 2 in LI2000. Each telescope is 
assumed to have a 10% total efficiency, and the data points 
are subject to photon counting errors only. The integration 
times for both telescopes are 100, 200 and 400secs for model 
(a), (b) and (c) respectively. 

The 10-m class telescope observations represent the cur- 
rent astronomical technology. The fractional deviation in- 
troduced into a fight curve by model (a), of order ~ 2%, is 
clearly detectable, with the photometry accurately reveal- 
ing the differing fight curve profiles at tfie different plane- 
tary orientations. This is also seen with model (b) which, 
while possessing lower signal-to-noise, still allows recovery 
of tfie detailed ligfit curve profile. For model (c), fiowever, 
the small fraction of reflected figfit means tfiat tfie resultant 
fractional deviations are small, less tfian 0.5%. As illustrated 
in Figure |^ wfiile fluctuations at tfiis level are potentially 
detectable, tfie details of tfieir profiles are lost in tfie noise. 
Plans are currently underway to build tfie next generation of 
ground-based telescope, witfi 30-m telescopes entering oper- 
ation between 2010-2015. An examination of tfie lower pan- 
els in Figure |^ demonstrate tfiat witfi tfie advent of sucfi 
telescopes, tfie fractional deviation introduced by a planet 
planetary companion will be readily identifiable, even for 
planets presenting only a tfiin crescent. 

It must be remembered tfiat tfie above discussion is for 
a particular set of planetary and microlensing parameters, 
sucfi as planetary radius, albedo, caustic strengtfi and rel- 
ative velocities, and tfiat tfie detectability of a planetary 
companion will depend on tfie specific values of tfiese for 
the system under consideration. 



4 DISCUSSIONS & CONCLUSIONS 

Recent studies fiave demonstrated tfiat wfien a star in tfie 
Galactic Bulge or Magellanic Clouds is magnified during 
a microlensing event, any figfit reflected from an orbiting 
planet can also be signiflcantly magnified to observable lev- 
els. Tfiese previous works, fiowever, fiave assumed a simpli- 
fied model for tfie distribution of reflected figfit from tfie 
planet, namely a circularly symmetric disk. The refiected 
figfit from any real planet, fiowever, will display phases akin 
to those seen witfi tfie moon and tfie inferior planets. 

Tfiis paper fias presented a study on tfie effect of plane- 
tary pfiase on tfie degree to wfiicfi tfie ligfit reflected from a 
planet is magnified as it is swept by a caustic during a grav- 
itational microlensing event. Considering simple models to 
represent various planetary pfiases, tfiis work fias demon- 
strated tfiat tfiese non-circular profiles result in strikingly 
different ligfit curves tfiat are solely dependent upon tfieir 
orientation witfi respect to tfie sweeping caustic. Wfiile tfiis 
befiaviour is different, depending upon tfie specifics of tfie 
ligfit distribution reflected by tfie planet, all tfie models dis- 
play tfie same generic trends. 

Equally interesting is tfiat tfie peak magnification of a 
ligfit curve for a model of tfie planet is also strongly de- 
pendent upon the planetary orientation witfi respect to tfie 
caustic. In eacfi model, tfie planet was substantially more 
magnified if orientated at ±7r/2, tfian at (see Figure 
In the case of the simple 'quarter moon' (model a), tfiis re- 
sulted in up to ~ 60 per cent more magnification tfian wfien 
compared to planets orientated at ~ 0. Due to scaling re- 
lations. Figures ^ and ^ are insensitive to specific values of 
caustic strengtfi and planetary radius. 

Tfiese trends are also seen in models (b) and (c) wfiicfi 
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Figure 6. A comparison of the flux received during the microlensing of a particular planet at differing phases and orientations. The 
orientation of the planet with respect to the incoming caustic follows the convention of Figure |l| and the phases under consideration 
are presented in the inset box. All the light curves are normalized to the peak flux received when the microlensed planet is full (i.e. a 
circular source), which corresponds to 100% on these plots. All panels display the same feature, that the amount of flux received falls 
monotonically as we consider more crescent- like source (left to right in the inset box), even though at these phases the planet is more 
strongly magnified. 



are more crescent-like, both possessing peak magnification 
distributions that possess a minimum at S ~ 0, rising to 
±7r/2. In the case of model (c), the light curve at —it/2 
possesses a peak magnification exceeding those at ^ ~ 
by 140 per cent. Considering that model (a) at 6' = is 
degenerate with the circularly symmetric models typically 
employed in microlensing studies, these results imply that 
previous results have underestimated the potential magni- 
fication in planetary microlensing events by a factor of a 
few. 

While these results are encouraging, the detectability of 
a planetary microlensing event will depend on the amount 
of flux received by an observer. As well as the magnification, 
this will depend also on the amount of stellar light that is re- 
flected by a planet, a value that decreases as a planet appears 
more crescent-like. Simulated observations presented in this 
paper reveal that the microlensing of moderately crecsent- 
like planets (model b) are detectable with current telescopes 
at a temporal sampling that allows the detailed recovery of 
the underlying proflle. With more extreme crescents, the 
planetary event the fractional deviation of ^ 0.5% occurs 
too quickly to be effectively sampled by photometric moni- 
toring with a 10-m class telescope, although proposed 30-m 
class telescopes will resolve the fluctuations. The exact de- 
tails of the microlensing light curve deviations due to the 
presence of a planetary companion depend upon a number 
of parameters, although the results presented here and in 



previous studies (Graff & Gaudi 2000; LI2000) suggest that 
this approach should successfully uncover hot Jupiters at 
kiloparsec-scale distances. 

Due to condensed particulate matter in the atmospheres 
of planets, the light they reflect is expected to be polarized 
(Seager et al. 2000) , although the degree of this polarization 
is extremely small (~ 10~^). LI2000 considered the effect 
of microlensing magniflcation on this polarization signature 
and found that it too can be boosted to detectable levels, 
although at a fraction of a percent, the measurement of such 
polarizations still presents observational difficulties. In light 
of the results of the study presented in this paper, we have 
begun an extensive investigation into the further effects of 
the microlensing of planets, especially the expected polar- 
ization signatures from crescent-like sources. 
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Figure 7. Simulated observations for the light curves presented in Figures 2 & 4. For each model, presented vertically and labelled 
(a), (b) & (c), there arc two sets of panels, the upper panel representing observations with a 10-m telescope, while the lower panel are 
observations with a 30-m telescope. 
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